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ABSTRACT 

New theoretical evaluations of the RGB luminosity function 'bump' and the ZAHB luminos- 
ity covering the range of metallicities typical of galactic globular cluster are presented. The 
variation of the theoretical RGB bump and ZAHB levels due to the metallicity, original helium 
content, mixing length value, age, mass loss, bolometric corrections, opacities and equation 
of state adopted in the evolutionary models is also discussed. These new prescriptions have 
been taken into account for casting light on a longstanding astrophysical problem connected 
with the Red Giant Branch evolutionary phase, namely the discrepancy between the obser- 
vational and the theoretical luminosity of RGB bump. A sample of globular clusters with 
accurate evaluations of the bump luminosity and spectroscopical metallicity determinations 
has been selected. The Zero Age Horizontal Branch luminosity at the RR-Lyrae instability 
strip has been evaluated as accurately as possible, and the observational luminosity difference 
between the RGB bump and the ZAHB has been compared with the theoretical values. It 
is shown that there is no significant disagreement between observations and canonical stellar 
models. The possible applications of this result are also briefly discussed. 
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'Z, 1 INTRODUCTION 

X 

5_j The existence of a well developed and populated Red 

C^ Giant Branch (RGB) is one of the main characteristics of 
the Color Magnitude Diagrams (CMDs) of galactic Glob- 
ular Clusters (GCs). While the Teff location of the RGB 
provides a means to calibrate the mixing length parameter 
for GC stars (ChiefB, Straniero & Salaris 1995, Salaris & 
Cassisi 1996 hereinafter Paper I), the Luminosity Function 
(LF) of the RGB of GCs is a useful tool to test the in- 
ner chemical composition stratification, since the hydrogen 
abundance outside the degenerate helium core sampled by 
the thin hydrogen burning shell, affects the rate of evolution 
during the RGB phase. In particular, due to the fact that 
the shell is extremely narrow in mass (thickness of the order 
of 0.001-0.0001 Mq), its crossing any discontinuity in the 
hydrogen profile leads to a temporary drop in the luminos- 
ity and therefore to a characteristic 'bump' in the observed 
differential LF. 

It has been known for a long time (see e.g. Thomas 
1967, Iben 1968, Renzini & Fusi Pecci 1988, Castellani, Chi- 
effi & Norci 1989, Bono & Castellani 1992) that theoretical 
RGB LFs display a characteristic bump due to the passage 
of the thin hydrogen burning shell through the composition 
discontinuity left by the deepest penetration of the convec- 



tive envelope. In spite of this, the first clear identification 
of the RGB bump in the observed LF of GCs is relatively 
recent; it goes back to the work by King et al. (1985) on 
47 Tuc, and only more recently Fusi Pecci et al. (1990, 
hereinafter FP90 - 11 clusters), Bergbush (1993 - NGC288), 
Sarajedini & Norris (1994 - 47 Tuc plus other 5 clusters), 
Sarajedini & Forrester (1995 - NGC6584), Brocato et al. 
(1996 - NGC5286, NGC6266, NGC6934, NGC6981) have 
extended the sample of GCs with a detected RGB bump. 

The reason for the difficulty in detecting the bump from 
GCs photometry is the need of very large RGB star samples 
in order to obtain a firm identification from the observed 
LF. According to the discussion in FP90 at least 1000-1500 
stars in the upper 4 magnitudes of the LF are required for 
a clear identification of the bump in metal poor GCs; this 
star sample is larger than the ones currently available for the 
best studied metal poor clusters. At higher metallicities the 
extension in luminosity of the bump is larger, and it is also 
shifted to lower luminosity, in more populated RGB regions, 
so it is easier to detect. 

In order to avoid uncertainties due to the calibration of 
the observational data and to the cluster distance modulus 
it is more reliable to consider, when comparing theoretical 
RGB bump luminosities with the observations, the parame- 
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ter AV|j;j™P = (Vbump - Vhb) (FP90) defined as the differ- 
ence in visual magnitude between the RGB bump and the 
Zero Age Horizontal Branch (ZAHB) level. 

FP90 tried to use the bump luminosity as a standard 
candle, in order to calibrate the relation between luminos- 
ity of the HB and metallicity. Since the bump luminosity 
depends also on the age of the stellar populations, and both 
bump and ZAHB are affected by the He abundance, they 
explored different scenarios about the age of the GC popu- 
lation and the He content. They compared the observed val- 
ues of AV, JJ™'' from their sample with the theoretical ones 
obtained from RGB models by Rood & Crocker (1989) and 
unpublished HB models by Rood (transformed to the obser- 
vational plane by adopting the Bell & Gustafsson 1978 and 
Kurucz 1979 transformations); they derived that the best 
agreement was obtained by assuming a constant age (t=15 
Gyr) and a constant He abundance (Y=0.23), or a constant 
age (t=15 Gyr) and an He abundance scaling as dY/dZ — 3 
(according to Peimbert & Torres- Peimbert 1977). In these 
two cases the results of the comparison were: 
i) the run of the theoretical relation with respect to the 
metallicity is in very good agreement with the observations, 
ii) the absolute value of the theoretical relation is 0.4 mag 
too bright. 

A possible interpretation of this discrepancy has been pre- 
sented by Alongi et al. (1991); they interpreted the detected 
difference in the zero-point as a limitation of the standard 
stellar models, and proposed undershooting from the base of 
the formal boundary of the convective envelope to reconcile 
theory and observations. More recently Straniero, Chieffi 
& Safaris (1992) and Ferraro (1992) pointed out that the 
proper inclusion of the a-element enhancement in the origi- 
nal chemical composition of the GC stars (see e.g. Wheeler 
et al. 1989) could help in reducing only partially the dis- 
crepancy of the zero-point without invoking undershooting, 
but until now the question remains unsettled. 

In this paper we present new theoretical determinations 
of the RGB bump and ZAHB luminosities covering the range 
of metallicities typical of galactic GCs, and discuss the in- 
fluence of different physical and chemical inputs adopted in 
the model computations on the value of AV, u™''. Moreover, 

we reexamine critically the observational values of AV. y™^, 
and by comparing theory with available observational data, 
we will show that actually a significant discrepancy between 
standard stellar models and observations does not exist. 

The plan of the paper is as follows: the theoretical mod- 
els are presented in section 2, while section 3 deals with the 
discussion of the observational data and their comparison 
with the theoretical values of AV. y™^; conclusions follow 
in section 4. 



2 THEORETICAL DETERMINATION OF THE 
HB AND RGB BUMP LUMINOSITY LEVEL 

In order to compare the observational values of AV. u™^ 
with the theoretical prescriptions, we have computed canon- 
ical evolutionary models (no diffusion, no undershooting 
from the base of the convective envelope) of stars with 
masses of O.75M0, O.SOMq and O.9OM0 - covering a range 
of masses presently evolving along the RGB of GCs corre- 
sponding to ages between approximately 22 Gyr and 12 Gyr 



- metallicities Z^O.OOOl, 0.0003, 0.0006, 0.001, 0.003, 0.006, 
and Y=0.23. We have interpolated between these models 
for each metallicity, in order to obtain the value of Vbump 
corresponding to an average age t=15 Gyr for the clusters 
(the choice of this average age comes out from the extensive 
analysis about GCs ages by Chaboyer, Sarajedini & Demar- 
que 1992, and Salaris, Chieffi & Straniero 1993). Since until 
now it is not clear if an age spread really exists in the bulk 
of the GCs population, we prefer to adopt this conservative 
hypothesis about the age of the clusters, at least for the 
purposes of this paper (but see also section 4). This corre- 
sponds to our reference scenario, in which all the GCs are 
coeval and with the same He content. 

The He core mass and the surface He abundance at 
the He ignition corresponding to the same age of 15 Gyr 
have been adopted in order to compute Zero Age Horizontal 
Branch (ZAHB) models of different total masses for each 
value of Z. As for the calibration of the superadiabatic en- 
velope convection, the mixing length calibration described 
in Paper I, obtained by fitting the empirical Tgfj values 
of GC RGBs obtained by Frogel et al. (1983), has been 
adopted. The evolutionary tracks and the ZAHB models 
have been translated to the observational plane by adopting 
the Kurucz (1992) transformations. The theoretical values 
of AV. u '"'^ have then been computed by considering the dif- 
ference in V magnitudes between the mean V magnitude of 
the RGB bump region, and V^ahb taken at log(Teff ) = 3.85, 
corresponding approximatively to the average temperature 
in the RR Lyrae instability strip. 

The theoretical relation AV. JJ™^ - [M/H], obtained by 
means of the models previously described, is our reference 
relation that will be compared with the observations. In 
Table 1, the value of Vbumpi ^zahb ^^ log(Toff) = 3.85 and 
AV. JJ™^ corresponding to our reference case, are listed for 
the different metallicities considered. 

All the theoretical models have been computed adopt- 
ing the FRANEC evolutionary code (see Chieffi & Straniero 
1989). The OPAL opacity tables (Rogers & Iglesias 1992, 
Iglesias, Rogers & Wilson 1992) for T > lOOOOiC and the 
Alexander & Ferguson (1994) opacities for T < lOOOOiC have 
been used. Both high and low temperature opacity tables 
are computed adopting the solar heavy elements distribu- 
tion (Grevesse 1991). The electronic conduction is treated 
according to Itoh et al. (1983). The equation of state 
(EOS) by Straniero (1988) has been used, supplemented by 
a Saha EOS at lower temperatures, as described by Chieffi 
& Straniero (1989). 

Before starting the comparison between our theoretical 
determinations for AV. JJ™^ and the observations, we have 
also analyzed in detail (computing additional stellar models) 
the dependence of this quantity on other different physical 
and chemical inputs adopted in computing stellar models 
and on the transformations used to transfer the models from 
the theoretical plane to the observational one. Whereas in 
literature it is easy to find a lot of exhaustive investigation 
on the influence on the luminosity level of the HB due to the 
various inputs adopted in evolutionary computations (Caloi, 
Castellani & Tornambe 1978, Dorman, Rood & O'Connell 
1993 and references therein), a detailed investigation con- 
cerning the properties of the RGB bump and therefore of 
the AV. u™'' is, till now, still lacking. In the following we 
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8) bolometric corrections. 



2.1 Metallicity and heavy elements distribu- 
tion. 

The luminosity of the RGB bump is strongly affected 
by a variation of the global amount of heavy elements. As a 
matter of fact a metallicity increase causes a decrease of the 
luminosity level of the bump. This occurrence is due to the 
larger extension in mass of the convective envelope during 
the first dredge up, related to the larger opacity and then, 
larger values of the radiative gradient with respect to the 
adiabatic one. Therefore the H discontinuity, located deeper 
and deeper with increasing metallicity, is reached by the 
H shell burning earlier during the RGB evolution. Fitting 
the values of Vbump in Tab. 1 as a function of [M/H] = 
log(M/H)star - log(M/H)0 fa log(Z) + 1.65 (where M here is 
the global metal abundance and Z the global heavy elements 
fraction) for t=15 Gyr yields: 



^bump^ 2.212 -h 1.768' 



[M/H] + 0.294 • [M/H]^ (1) 



with a r.m.s.=0.02 mag. 

Also the luminosity of the ZAHB is strongly affected 
by the heavy elements abundance, due to the variation of 
the He core mass at the flash with the metallicity. The 
higher the metallicity, the lower the He core mass and the 
ZAHB luminosity. From our reference models we derived 
the following relation between My'^'' at log(Teff) = 3.85 
and [M/H]: 



^zahb ^ 1.129 + 0.388- 



[M/H] + 0.063 • [M/H]^ (2) 



Figure 1. Top: Absolute luminosity of field and globular clus- 
ter RR-Lyrac stars versus metallicity provided by Clementini et 
al. (1995). The luminosity level of the ZAHB at the average 
temperature of the instability strip for two different assumptions 
concerning the equation of state are also displayed. The solid line 
was evaluated taking into account the EOS provided by Straniero 
(1988), whereas the dotted line is based on the EOS recently pro- 
vided by the OPAL group; Bottom: the same as the top panel 
but the absolute luminosities of the RR Lyrac variables have been 
estimated by adopting a different value for the conversion factor 
p (i.e. the parameter used for trasforming the radial velocities 
into pulsational velocities). 

will discuss the variation of AV. JJ™^ induced by changes of 
the following parameters: 

1) the metallicity and the distribution of the heavy elements; 

2) the abundance of Helium Y; 

3) the mixing length parameter; 

4) the age of the stellar population; 

5) mass loss during the RGB evolution; 

6) the equation of state; 

7) the opacity; 



with a r.m.s. =0.011 mag. 

In Fig.l a to b, we compare this relation with the most 
recent determination of the relation between absolute lu- 
minosity of field and globular cluster RR-Lyrae stars and 
metallicity, published by Clementini et al. (1995), based 
on new spectroscopical determinations of metallicity and on 
Baade Wesselink estimates of the absolute magnitudes of 
the variables. It is worth noting that in Fig.l, only observa- 
tional data for supposedly unevolved RR Lyrae stars have 
been considered and therefore they should provide a good 
estimation of the ZAHB luminosity at the RR Lyrae insta- 
bility strip. Moreover, Clementini et al. (1995) provide two 
different tabulations (their Tab. 21, columns 5 to 6) concern- 
ing the absolute luminosity of their RR Lyrae stars sample, 
corresponding to two different prescriptions for the conver- 
sion factor p between observed and true pulsational velocity. 
In panels a to b we have displayed their data corresponding 
respectively to columns 6 and 5 of their Tab. 21. We have 
derived the global metallicity [M/H] for the stars in their 
sample by assuming the [Fe/H] and [a/Fe] values given in 
the paper, and also the errors on the observational determi- 
nations of metallicity and luminosity come from their paper. 
As is evident from the figure, our theoretical relation is in 
agreement - for both choices concerning the conversion fac- 
tor p - with the observational data. 

From relations (1) and (2) we derive: 

AVJ'^P = 1.083 + 1.380 ■ [M/H] + 0.231 • [M/H]^ (3) 
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with a r.m.s. =0.023 mag. 

These relations, and in particular the third one, that 
we will use for the comparison with the observations, have 
been derived from stellar models computed adopting a scaled 
solar heavy elements distribution. As well known the orig- 
inal chemical composition of GCs stars is characterized by 
[a/Fe] > (see e.g. Wheeler et al. 1989); Salaris, ChiefR 
& Straniero (1993) have demonstrated that the evolution of 
low mass low metallicity stars with an a-enhanced heavy el- 
ements distribution and a fixed global metallicity [M/H], is 
very well reproduced by scaled solar models with the same 
value of [M/H]. 

We have verified that the same holds if we adopt the re- 
cent a-enhanced OPAL and Alexander & Ferguson molec- 
ular opacities (see Salaris et al. 1996). In particular the 
value of AV. JJ™^, obtained using these opacity tables for a 
fixed [M/H], is coincident - within 0.01 mag - with the one 
derived from the scaled solar models with the same [M/H]. 



2.2 The Helium abundance. 

The infiuence of a variation of the original He content on 
the value of AV, y™'' has been also tested. As well known, 
increasing the original He content increases the ZAHB lu- 
minosity, due to the more efficient energy generation in the 
H-burning shell. At the same time, also the level of the 
bump is shifted to higher luminosity, and the net effect is 
a slight reduction of the AV|jJJ™P. For each AY = -1-0.01 
we derive from our models a reduction of 0.01 Imag of the 
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2.3 The mixing length parameter. 

The mixing length parameter ml is one of the free pa- 
rameter which enters in stellar computations. However, as 
already discussed in Paper I, it is usually constrained by the 
requirement that the solar radius and the observed colors of 
red giant stars have to be correctly reproduced. In Paper 
I we have demonstrated that the value of ml obtained by 
reproducing the T^g of the Sun may not be suitable also 
for reproducing the T^g of the GCs RGB. For this reason 
we have performed the ml calibration for our evolutionary 
tracks of metal poor low mass stars by reproducing the ob- 
servational Toff of GCs RGBs as derived by Frogel et al. 
(1983). 

In order to illustrate the dependence of the bump lumi- 
nosity on the value of the mixing length, we have reexamined 
the evolutionary tracks discussed in Paper I, computed for 



different values of ml. We find that 



AVk 



Ami 



'0.27mag. 



Since we have verified that the level of the ZAHB in the RR 
Lyrae stars region is not infiuenced (as it is well known) by 
a variation of ml, it is obtained that the value of AV, y™'' 
is decreased by ~ 0.04 mag for a variation by -1-0.15 of ml. 
It is worth noting that a variation of ml of about 0.15 
is large enough to compromise the agreement between theo- 
retical evolutionary models and theT^ff of GCs RGB stars. 
In Paper I we derived a relation between [M/H] and Toff of 
the GCs RGB by adopting the Tpff determinations by Fro- 
gel et al. (1983); this relation has been used for calibrating 
all the stellar models presented in this paper. The disper- 
sion of the observational points around this relation is of the 



order of lOOK (see the discussion in Paper I) , and since this 
dispersion corresponds to a variation of the ml by around 
±0.10, we can assume this quantity as an estimate of the 
maximum uncertainty associated to the calibrated value of 
ml. Taking into account this indetermination, the variation 
of the theoretical values of AV. J^™^ results to be less than 
±0.03 mag, that is quite negligible. 



2.4 The age of the stellar population. 

The location of the bump on the RGB of a stellar track 
depends on the mass of the model, i.e. on the age of the stel- 
lar system in which that star is now evolving (see Straniero 
& Chieffi 1991), while the ZAHB luminosity is practically in- 
dependent on the age, at least in the age range spanned by 
the GCs. From the evolutionary tracks of O.75M0 - O.8M0 
- O.9M0 models, we derive for each metallicity an increase 
of AV. 1^™'' by almost 0.024 mag for an increase of 1 Gyr in 
the age of the clusters. 



2.5 The effect of the mass loss on the RGB. 

As it is well known, all the low mass stars during their 
evolution along the RGB experience the phenomenon the 
mass loss. Mass loss during the RGB is also required to ob- 
tain the correct HB morphology of GCs. However mass loss 
does not affect at all the main evolutionary properties of the 
stars, as for instance the helium core mass at the he flash 
and the location in Toff of the RGB. This occurrence is re- 
ally correct if "conventional" assumptions are made on the 
efficiency of this phenomenon. In fact, it has been shown 
(Castellani & Castellani 1991, D' Cruz et al. 1995) that, 
assuming a very high efficiency of the mass loss mechanism, 
it could be possible to "obtain" stellar models whose evolu- 
tion does not follow the prescriptions of the canonical stellar 
evolution theory. 

We have tested if mass loss could affect in some way 
the bump luminosity on the RGB. For this aim, some evo- 
lutionary tracks have been computed adopting various as- 
sumptions on the mass loss efficiency. Let us remember that 
the mass loss phenomenon is usually parametrized in stellar 
computations the Reimers (1975) formula, in which appears 
a free parameter 77. According to various authors (see for 
instance Renzini & Fusi Pecci 1988) to finely reproduce the 
distribution of stars on the HB of the bulk of galactic glob- 
ular cluster - except the GCs showing HB blue tails in their 
CMDs it is necessary to adopt for such a parameter a value 
around 0.3 - 0.4. The evolutionary tracks corresponding to 
the same O.8M0 model computed without mass loss in a 
case and assuming a large efficiency (r; = 1.0) for the mass 
loss are plotted in Fig. 2. One can easily notice that the 
effect of the mass loss on the bump is absolutely negligible. 
It is also worth noticing that this result has been obtained 
assuming a very strong amount of mass loss, as it is con- 
firmed by the occurrence that the star is forced to leave off 
the RGB before igniting the He central burning. Therefore 
one can safely assume that, in the range of efficiency of the 
mass loss phenomenon in real RGB stars, the influence of 
this mechanism on the bump luminosity is quite negligible. 

2.6 The equation of state. 
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GCs. 

Once calibrated the mixing length parameter on the ob- 
servational data by Frogel et al. (1983), we have computed a 
set of evolutionary tracks as described for our reference case; 
ZAHB luminosities result to be about 0.055 mag higher than 
for the models computed with the Straniero EOS (see Fig.l) 
and AV. y™^ values are 0.04-0.05 mag lower than the ones 
obtained adopting our reference scenario and the Straniero 
EOS. When adopting an age of 15 Gyr also for the models 
computed with the OPAL EOS, AVJ^™^ values 0.03 mag 
higher than the reference values are obtained. 



Figure 2. The H-R diagram for a O.SMq stellar model computed 
under two different assumptions concerning the efficiency of the 
mass loss phenomenon. In the inset the effect of mass loss on the 
evolutionary tracks in the RGB bump region is shown. 

A new equation of state suitable for stellar evolutionary 
computations has been very recently provided by Rogers, 
Swenson & Iglesias (1996) (OPAL EOS). Therefore we have 
decided to test if theoretical values of AV. JJ™^ are modified 
when using this new physical input in the evolutinary models 
computations. In this case we have supplemented the OPAL 
EOS (in the regions not covered by the tables), with a Saha 
EOS (for T < ^QQQK) and the Straniero (1988) EOS, as 
described in Salaris, DeglTnnocenti & Weiss (1996). 

The use of the OPAL EOS in evolutionary computa- 
tions has led to a revision of the age of the GCs with respect 
to determinations obtained by means of stellar models com- 
puted using a different EOS (as in Chaboyer et al. 1992 
and Salaris et al. 1993). Recent works by Chaboyer & Kim 
(1995, who find a reduction by 6-7% in the ages derived us- 
ing My {TO) and an average age of 13 Gyr for a sample of 
40 clusters), MazziteUi, D'Antona & Caloi (1995), Salaris et 
al. (1996), show that the average age of the GCs should be 
around 12-13 Gyr, significantly lower than the average age 
obtained in previous works. This result is very important 
since, as it is well known, the determination of the GCs age 
is a fundamental tool to investigate the galactic formation 
mechanism and the age of the Universe. The discussion of 
these results is out of the goals of the present work, but 
due to the effect of the age on the RGB bump location, 
we have to make some realistic assumptions concerning the 
adopted value for the age of the GCs. Using the same cri- 
teria adopted to obtain our reference models, in computing 
the theoretical values of AV, J^™'' by using the OPAL EOS 



2.7 The opacity. 

The bump luminosity, being related to the position of 
the H discontinuity produced by the convective envelope 
during its deeper penetration, depends strongly on the opac- 
ity evaluation for temperatures of the order of one milion of 
degrees, i.e. the temperature at the bottom of the convective 
envelope. This occurrence has been tested in the past when 
changing the opacity tables from the Cox & Stewart (1970) 
and Cox & Tabor (1976) to the most reliable Los Alamos 
opacity library (hereinafter LAOL; Huebner et al. 1977) the 
brightness of the bump decreased by ~ 0.2mag; the change 
from the LAOL to the OPAL opacities causes a much smaller 
reduction of the bump luminosity, by ~ 0.07mag. Now the 
opacity evaluations are more accurate than in the past, both 
in the high temperature and in the low temperature region 
(see e.g. Paper I for a comparison between three recent sets 
of low temperature opacities). So one can be hopefully confi- 
dent that there is not much room for a significative variation 
in the bump brightness as due to a forthcoming generation 
of updated opacity libraries. For instance, in very recent 
time, a big effort has been made to improve the accuracy of 
the OPAL opacity evaluations, by increasing the number of 
elements taken into account in the metal mixture (Roger & 
Iglesias 1995 - 21 elements mixture). The effect of this last 
generation of the OPAL opacities on the bump luminosity 
has been also checked. As a result, there are not significative 
variations of AV, !J™^. 

2.8 Bolometric corrections. 

Our theoretical evolutionary tracks have been trans- 
formed into the observational plane by adopting the Kurucz 
(1992) transformations. In order to check the sensitivity of 
AV, JJ™'' to different sets of bolometric corrections adopted, 
we have also used the transformations by Buser & Kurucz 
(1992) supplemented with the Buser & Kurucz (1978) ones 
for T > 6000/i", as described in Salaris et al. (1996). The 
same values of AV. u™'' - within 0.01 mag - have been ob- 
tained at each metallicity. 



3 THEORY VERSUS OBSERVATIONS. 

As discussed in the first section, RGB bumps have been 
detected only in a few globular clusters. To perform a 
meaningful comparison between theory and observations, we 
have to consider clusters with a quite clear detection of the 
bump, and with an accurate determination of the metal- 
licity (taking into account also the enhancement of the a- 

bump 



we have therefore chosen an average age of 12 Gyr for the elements), since the AV,, is strongly dependent on the 
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adopted global metallicity. For this reason only clusters with 
high resolution spectroscopical determinations of the photo- 
spheric abundances of Fe and a elements have been taken 
into account. In Tab. 2 the seven clusters considered, the 
adopted values of [M/H] (these values come from Safaris & 
Cassisi f996, where determinations of [a/Fe] and [Fe/H] are 
collected for a sample of 22 globular clusters) , the values of 
Vzahbi Vbump) '^nd AV. JJ™'' (with its associated observa- 
tional error) are reported. The sources of the photometric 
data are given in the discussion of the individual clusters 
(see below). Five clusters of this GCs sample are very poorly 
populated or not populated at all in the RR Lyrae stars re- 
gion (NGC6397, NGC6752, M79, 47 Tuc and NGC6352), 
and the determination of the ZAHB luminosity in the in- 
stability strip is not straightforward. The values of V^ahb 
reported in Tab. 2 have been carefully determined by adopt- 
ing different procedures, depending on the morphology of 
the cluster HB. 

In the case of M3 and M5, whose CMDs display a very 
well populated HB in the RR Lyrae region, we have adopted 
the mean luminosity of the RR Lyrae variables ((Vrr)) as 
provided in the original works (see below). Clearly (Vrr) 
does not represent the value of the ZAHB luminosity, that 
has to correspond to the lower envelope of the observed HB 
stars distribution (see e.g. Sandage 1990); thus, for obtain- 
ing the ZAHB level one has to correct the (Vrr) value tak- 
ing into account the thickness of the observed HB. Carney 
, Storm & Jones (1992) provide a relation between V^ahb 
and (Vrr) derived by using the data published by Sandage 
(1990), who carefully studied the vertical distribution of HB 
stars in various GCs. They used a sample of 8 clusters to 
derive the following relation: 



V. 



zahb 



(Vrr) + 0.05[Fe/H] + 0.20 (4) 



that gives the ZAHB luminosity as a function of the mean 
luminosity of the RR Lyrae stars and of the cluster iron 
content. We have performed the same kind of analysis by 
using clusters, in the sample considered by Carney et al. 
(1992), for which spectroscopical determinations of [Fe/H] 
and [a/Fe] (see Paper I) are available, thus obtaining the 
following relation: 

VzAHB = (Vrr) + 0.04[M/H] + 0.15 (5) 

Relation (5) has been used in order to correct the val- 
ues of (Vrr) for obtaining the ZAHB level. In each case 
this formula has been used, we have verified that the result- 
ing luminosity matches the lower envelope of the HB stellar 
population (see Fig. 3a,b). 

In the case of NGC6397, NGC6752 and M79 which are 
characterized by blue HBs, we have followed the same pro- 
cedure used by other authors (see Buonanno et al. 1986, 
Alcaino et al. 1987, Ferraro et al. 1992) in order to de- 
rive their ZAHB luminosity. We have considered a reference 
cluster with the same metallicity but with an HB populated 
in the blue and in the RR Lyrae region, and its CMD has 
been shifted in such a way that the blue part of its HB se- 
quence overlaps the blue HB of the cluster we are studying, 
in order to form a unique sequence (see Fig. 3c-e). In this 
way the ZAHB level of the cluster with the blue HB has been 
obtained from the reference cluster, after the correction for 



the relative luminosity shift. When it has been selected as 
reference cluster a GC with the appropriate metallicity, and 
if the CMD shift has been correctly done, the two RGBs 
have also to be perfectly overlapped (see Fig. 3c-e), and the 
horizontal shift applied to the reference cluster has to corre- 
spond to the reddening difference between the two clusters. 

As for 47 Tuc and NGC6352 a procedure very similar 
to that described in FuUton et al. (1995) has been adopted. 
The HBs of these two clusters are populated only on the 
red side of the instability strip, so no estimate of the ZAHB 
level at the region of the RR-Lyrae stars is available. More- 
over, the value of the ZAHB at the red side can not be used 
as an estimate of the ZAHB luminosity for the RR Lyrae 
region (see the discussion in Castellani, Chieffl & Pulone 
1991). In this case we have derived from our ZAHB models 
(transformed to the observational plane by adopting both 
the Kurucz 1992 and the Buser & Kurucz 1978, 1992 trans- 
formations) for Z=0.003 and Z=0.006 the difference 5 in My 
between the red part of the ZAHB and the point along the 
ZAHB at log(Teff ) = 3.85. We obtained S = 0.10±0.05, and 
applied this correction to the observational data, in order to 
perform the comparison with the theoretical luminosities at 
log(Teff ) = 3.85. 

In the following the case of each cluster (in order of 
increasing metallicity) will be discussed separately: 

(i) NGC6397: The luminosity of the bump (Vbump = 
12.60 ± 0.10) is derived from FP90 who adopt the photom- 
etry by Alcaino et al. (1987). In the paper by Alcaino 
et al. (1987) the mean luminosity of the HB in the RR- 
Lyrae stars region was obtained by superimposing the CMD 
of M15 (which is populated in the RR Lyrae region) to 
that of NGC6397 (see previous discussion). We have fol- 
lowed the same procedure (see Fig. 3c) using the CMD of 
M68 (Walker 1994), which has a global metal abundance 
much more similar to that of NGC6397 according to the 
latest spectroscopical determinations ([M/H]=-1.78, see Pa- 
per I), and a well populated HB in the variable stars re- 
gion. Walker (1994) provides (Vrr) = 15.64 ± 0.01 for 
the RR Lyrae stars in M68, to which we have subtracted 
a quantity AV — —2.70, corresponding to the shift we have 
applied to the M68 diagram in order to match the blue 
HB of NGC6397. The horizontal shift (A(B - V)) that 
we applied to M68 is A(B — V)=-(-0.1, in good agreement 
with the reddening difference between these two clusters 
(E(B - V)m68 = 0.06 - 0.10 according to Walker 1994, and 
E(B — V)i\fGC6397 = 0.17 — 0.20 according to Alcaino et al. 
1987). Estimating an uncertainty of around 0.1 mag due to 
the procedure previously described, after using relation (5) 
a value Vzahb ~ 13.02 ± 0.10 is finally derived. 

(ii) NGC5272 (M3): The observational data come from 
Buonanno et al. (1994). In this paper the authors provide 
the value of Vbump = 15.40±0.05 and (Vrr) = 15.66±0.05. 
By applying the correction for the ZAHB we obtain V^-ahb ~ 
15.76 ±0.05. 

(iii) NGC6752: The bump level (Vbump = 13.65 ± 0.05) 
comes from FP90 and has been determined by using the 
photometry by Buonanno et al. (1986). In order to de- 
rive the ZAHB level we have superimposed the data of M3 
(which has a very well populated RR Lyrae region) to that of 
NGC6752 (see Fig. 3d ), since the two clusters have quite the 
same metallicity (see Tab. 2). The ZAHB luminosity of M3 
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Figure 3. CMDs of the clusters for which the ZAHB luminosity level has been obtained using our relation (5), or superimposing 
another cluster populated in the RR Lyrae region (see text). In each panel the adopted ZAHB luminosity level and the associated error 
is indicated. In panel c) the C-M diagram of M68 (filled triangles) is superimposed to that of NGC6397 (open circles); the same in panel 
d) and e), but with M3 (filled triangles) superimposed respectively to NGC6752 and M79 (open circles). 



has been corrected by a quantity AV — —1.90 corresponding 
to the vertical shift applied to its CMD in order to match the 
data of NGC6752. The horizontal shift A(B - V) = +0.05 
applied to M3 agrees very well with the reddening differ- 
ence between the two clusters (E(B — V)m3 = 0.00 — 0.03 
according to Buonanno et al. 1994, and E(B — V)ngC6752 = 
0.02 — 0.06 according to Penny & Dickens 1986). We derive 
Vzahb = 13.86 ± 0.11, after taking into account an error of 
about 0.1 mag due to the procedure of superimposing the 



two HB sequence, and the error associated to the estimate 
of the M3 ZAHB level (see above). 

It is very interesting to note the agreement between the 
cluster distance modulus obtained by adopting this value 
of Vzahb together with our relation (2), and the distance 
modulus determined very recently by Renzini et al. (1996) 
in a completely independent way. They use the observed 
White Dwarf (WD) cooling sequence of NGC6752 as a dis- 
tance indicator, and derive the distance modulus by fit- 
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Table 2. Luminosity of the bump, of the Horizontal braneh at 
log(T(,ff) = 3.85 and AV, , ^ for the clusters considered. 

rbump 



Name 



V, 



zahb 



Vb 



unip 



AV 



hb 



[M/H] 



NGC104 


14.20 


14.55 


0.35 ±0.18 


-0.70 


NGC1904 


16.36 


16.00 


-0.36 ±0.12 


-1.27 


NGC5272 


15.76 


15.40 


-0.36 ±0.07 


-1.31 


NGC5904 


15.15 


14.95 


-0.20 ±0.07 


-1.19 


NGC6352 


15.50 


15.86 


0.36 ±0.12 


-0.70 


NGC6397 


13.02 


12.60 


-0.42 ±0.14 


-1.70 


NGC6752 


13.86 


13.65 


-0.21 ±0.12 


-1.28 



ting the cluster WD sequence to an empirical cooling se- 
quence constructed using local WDs with well determined 
trigonometrical parallaxes. Following this procedure they 
derive (m — M)o = 13.05 with an overall uncertainty less 
than ±0.1 mag. By adopting our relation (2) together 
withVzahb = 13.86 ± 0.11, and Ay = 0.12 ± 0.06 (see 
Renzini et al. 1996, Penny & Dickens 1986), we obtain 
(m — M)o = 13.01 ± 0.13, in good agreement with the result 
by Renzini et al (1996). 

(iv) NGC1904 (M79); The observational data come from 
Ferraro ct al. (1992). We adopt Vbump = 16.00 ± 0.05 as 
estimated by the authors. From the observational data wc 
have determined the ZAHB level following the same proce- 
dure as described for NGC6397 and NGC6752. Also in this 
case we have shifted the data of M3 in order to superimpose 
the two HB sequences (see Fig.3e). The vertical shift applied 
to M3 is AV = +0.6, while no correction to its (B-V) values 
has been applied; this is in good agreement with the fact that 
the reddening of the two clusters (E(B-V)m79 = 0.00-0.02 
according to Ferraro et al. 1992) is coincident. The ZAHB 
luminosity of M79 results to be V^ahb = 16.36 ± 0.11. 
(v) NGC5904 (M5): The bump luminosity has been 
derived from Brocato, Castellani & Ripepi (1995, 1996) 
(Vbump ~ 14.95 ±0.05). From the same authors we have de- 
rived an average magnitude for the RR Lyrae stars (Vrr) — 
15.05 ± 0.05. By applying relation (5) V^ahb = 15.15 ± 0.05 
is obtained. 

(vi) NGC104 (47 Tuc) and NGC6352; The bump 
(Vbump = 14.55 ± 0.05) and the ZAHB luminosities 
(Vzahb = 14.10 ± 0.15) of 47 Tuc come from FP90, who 
adopted the photometry by King et al. (1985). As for 
NGC6352 we have used the data from Sarajedini & Nor- 
ris (1994) obtaining V^ahb = 15.40 ± 0.10 from the fit 
of the lower envelope of the HB stellar population, and 
Vbump = 15.86 ± 0.05 as quoted by the authors. FuUton 
et al. (1995) published recently another photometry for 
NGC6352, obtained in part from ground based observations 
and in part from HST observations, and they found a zero 
point difference in the V magnitudes between their work and 
the results by Sarajedini & Norris, for the 145 stars in com- 
mon. However, this zero point difference (around 0.16 mag) 
does not affect at all - as obvious - the differential quantity 



bump 
hb 



AV- 

By applying the correction 5 — 0.10 ± 0.05 due to their 
red HB (see the previous discussion), we obtain V^ahb — 
14.20 ± 0.17 for 47 Tuc and V^ahb = 15.50 ± 0.11 for 
NGC6352. 

In Fig. 4, the comparison between the observational val- 
ues of AV, y™^ for the seven clusters considered and the 



Figure 4. The values of AV, , ^ versus the global metallicity 
for all the clusters in our sample. Our reference relation (see 
text) (solid line) and the theoretical prescription obtained using 
in stellar computations the OPAL EOS (dashed line) are also 
plotted. 

theoretical prescriptions is reported. We have taken into ac- 
count the observational errors in the determination of such 
a quantity (see Tab. 2), and assumed an uncertainty of 
0.15dex on the spectroscopic estimates of [M/H] (see the 
discussions in Gratton, Quarta & Ortolani 1986, Gratton & 
Ortolani 1989 and Kraft, Sneden, Langer & Shetrone 1993). 
The theoretical relations corresponding to the models com- 
puted adopting the Straniero (1988) EOS and the OPAL 
EOS are displayed. These two relations correspond - as pre- 
viously discussed - to an age of 15 Gyr (Straniero EOS) 
and 12 Gyr (OPAL EOS), and they have been computed by 
adopting Y=0.23. 

The figure clearly demonstrates the overall agreement 
between theoretical standard stellar models (computed with 
the Straniero or the OPAL EOS) and observations, either 
for the run of AV, J^™^ with respect to [M/H] , either for the 
absolute values of this quantity; all the observational points 
are fitted within the observational error bars. The worst 
agreement is obtained for NGC6397, the most metal poor 
cluster in our sample. 

We may therefore conclude that there is no signiGcant 
discrepancy between observations and canonical stellar mod- 
els computed by adopting updated input physics, under the 
assumption of a constant age for the GCs, and a constant 
initial Helium abundance. 



4 DISCUSSION AND CONCLUSIONS. 

In the two previous sections we have presented new 
theoretical stellar models computed with updated input 
physics, with the aim of comparing the theoretical values 
of AV. u™'' with the observations. Before comparing this 
quantity with real clusters, we had to make some assump- 
tions about the age and the original He content of our 



Red Giant luminosity function 'Bump ' 9 



Figure 5. As in Fig. 4, but the tiieoretical relation AV. , - 

[M/H] are now displayed for our reference case (solid line), witii 
two different assumptions concerning the cluster age (13Gyr - 
short dashed line; 17 Gyr - long dashed line) and for the case 
of an Helium abundance scaling with the nietallicity (see text) 
(dotted line). 

models, since the AV. ^""^ is influenced both by the age 
(a change of the age affects the bump luminosity) and by 
the value of Y (a variation of the Helium content modi- 
fies the ZAHB and bump luminosities). In the conservative 
hypothesis of coeval GCs with the same original He con- 
tent (Y=0.23), a very good agreement between theory and 
observations can be derived from Fig. 4. This means that 
standard stellar models can actually reproduce the luminos- 
ity levels of ZAHB and bump, and their run with respect to 
the metallicity; it confirms also that standard stellar models 
reproduce accurately the Hydrogen profile in the interior of 
RGB GCs stars. The differences between our results and 
the conclusions by FP90 are due basically to three different 
reasons: 

i) the use of updated evolutionary models; 

ii) the adoption of new spectroscopical determinations of 

[a/Fe] and [Fe/H]; 

iii) new observational data. 

FP90 have used models by Rood computed with the old Cox 
& Stewart (1970) and Cox & Tabor (1976) opacities and, as 
discussed in their paper, it is not taken into account in the 
computation of their HB models the extra helium brought 
to the surface during the first dredge up. In addition, as 
yet discussed in the introduction, they used also a differ- 
ent set of bolometric corrections (Bell & Gustafsson 1978 
and Kurucz 1979). The net resulting difference between our 
theoretical values of AV. JJ™^ and the ones used by FP90 
(that arises from the differences on both the ZAHB and the 
Bump levels) ranges from 0.12 up to 0.22 mag (it is higher at 
higher metallicities). The remaining part of the discrepancy 
between theory and observations disappears when consider- 
ing individual high resolution spectroscopical determination 



of [M/H], including the overabundance of the a elements, 
and adopting (when they are available) more recent obser- 
vational data. 

As a further step we can now check how the fit pre- 
sented in Fig. 4 is modified by relaxing the hypothesis of co- 
eval GCs or of a constant Y. If we consider a constant Y but 
an age spread by ±2 Gyr around our assumed average age 
of the clusters (see e.g. Chaboyer, Demarque & Sarajedini 
1996), the theoretical AV. JJ™^ values would be spread by 
approximately ±0.05 mag around the line displayed in Fig.4. 
In Fig. 5 two lines corresponding to the AVj^JJ™^ values ob- 
tained for our average age increased and decreased by 2 Gyr 
are displayed. As for a variation of Y with the metallicity, 
the effect of a dY/dZ = 3 on the AV, v™^ is shown in Fig. 5. 
It is evident from the figure that the adopted variation of Y 
with the metallicity has a negligible effect on the theoreti- 
cal AV[J™P values (at most 0.02 mag at [M/H] =-0.6), and 
also the assumed spread in the ages of the clusters does not 
change very much the theoretical AV, y™^-[M/H] relation, 
remaining always compatible with the observational data. 
Because the agreement between our theoretical AV. JJ™^ val- 
ues and the observations, relation (3) can be safely used for 
deriving the clusters metallicities. By simply differentiating 
this relation, one obtains that a variation of the metallicity 
by ±0.15 dex around [M/H] = — 1.3 corresponds to a varia- 
tion of AV. u™'' by ±0.12 mag. Allowing for an observa- 
tional error of around 0.12 mag (see for example Tab. 2) in 
the determination of AV. u™'', and for an error of around 
0.1 mag due to uncertainties in the age, He content of the 
cluster and the mixing length calibration of the theoretical 
models (see section 2), it would be possible to estimate the 
metallicity of the cluster with an error of around ±0.20 dex 
(if [M/H] is around -1.3). It is worth noticing that the lu- 
minosity difference between the bump and the Horizontal 
branch is a metallicity index which relies only on the obser- 
vation of luminous stars in the GCs and is less dependent 
on the precise value of the mixing length than others metal- 
licity indicators based on the colors or on the shape of the 
RGB. 

When the metallicity of a cluster is known, it is possi- 
ble to use the bump luminosity as a standard candle. Re- 
lation (3) can be used as a tentative guess to estimate the 
HB luminosity level in those globular clusters in which the 
horizontal branch is poorly populated near the RR-Lyrae 
instability strip or the HB is quite blue or quite red. It 
could be interesting to test this procedure in evaluating the 
Helium abundance on the basis of the R method in those 
clusters for which estimating the HB luminosity level is a 
thorny problem. 

Before concluding we want to stress again that in this 
investigation we have computed canonical stellar evolution- 
ary models, neglecting other non canonical effects, as, for 
instance, the Helium (and heavy elements) diffusion. The 
inclusion of this mechanism in solar models changes, among 
other quantities, the chemical abundances in the convective 
envelope and the depth of the convective region with respect 
to canonical models, improving the agreement with helioseis- 
mological data (see, e.g., the discussion in Castellani et al. 
1996). 

Proffit & Vandenberg (1991) have studied the evolution 
of GCs stars taking into account He diffusion, and very re- 
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cently Castellani et al (1996) have shown that He and heavy 
elements diffusion does not change appreciably (age differ- 
ences by less than 1 Gyr) the age of the GCs with respect 
to canonical evaluations, but at present time an analysis of 
the influence of the He and heavy elements diffusion on the 
AV. u™^ does not yet exist; from preliminary computations 
at Z=0.0004 (Degl'Innocenti 1996, private communication) 
it results that this quantity is changed by only ~ 0.02 mag 
with respect to the canonical value, but in order to test 
further the stellar models which include diffusion, computa- 
tions also for higher metallicities have to be performed. 
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